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ABSTRACT 

We  define  the  door-width  of  a  simple  polygon  (a  'chair')  and  give  an 
O(n-)  algorithm  for  computing  its  door-width.  It  is  first  shown  that  all 
passages  of  the  chair  through  the  door  can  be  reduced  to  a  sequence  of 
certain  elementary  motions.  We  introduce  the  technique  of  constraint 
analysis  in  characterizing  elementary  motions.  Our  algorithm  actually 
constructs  a  motion  of  the  chair  through  a  door,  and  thus  is  a  'local- 
expert'  for  planning  motion  through  doors.  Such  algorithms  have 
applications  in  more  general  motion-planning  systems  in  robotics. 
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1.    Introduction 

To  move  a  given  chair  through  a  door,  a  natural  question  is  how  wide  the  door 
should  be.  The  minimum  width  of  a  passable  door  is  called  the  width  of  the  chair.  The 
question  was  first  raised  by  Gilbert  Strang  [Str82]  who  studied  the  case  of  convex  chairs 
in  two  and  three  dimensions.'  We  focus  on  the  polygonal  chairs  in  this  paper  but  allow 
non-convex  chairs.  Although  the  notion  of  width  for  convex  chairs  is  easily  captured,  the 
width  of  non-convex  chairs  admits  several  reasonable  interpretations.  We  will  note  three 
possibilities  in  this  paper. 

To  compute  the  width  of  a  chair,  we  in  fact  construct  a  motion  of  the  chair  through 
any  door  of  sufficient  width.  Thus  our  problem  can  be  regarded  as  an  instance  of  the 
motion-planning  problem  that  has  been  extensively  studied  in  the  robotics  literature.  We 
refer  to  [Yap85]  for  a  survey.  It  should  be  observed  that  our  problem  is  a  curious  one 
from  the  robotics  point  of  view:  in  robotics,  one  is  typically  interested  in  algorithms  for 
moving  some  fixed  robot  B  through  an  arbitrary  environment  E  (from  some  well-defined 
class  of  environments)  where  E  is  described  as  part  of  the  problem  instance.  In  our  case, 
we  have  an  algorithm  for  a  fixed  environment  E  but  permit  the  robot  B  to  be  arbitrarily 
specified  (from  the  class  of  polygonal  robots).  In  [Yap85]  we  term  the  former 
basic  planning  problems  and  the  latter  local  planning  problems .  A  motion-planning 
algorithm  that  is  specialized  for  a  fixed  locality  E  is  thus  a  "local  expert'.  Although  the 
papers  [HJW82],  [KK84]  were  not  originally  so  conceived,  we  can  regard  them  as  further 
examples  of  local  planning. 

We  indicate  the  pragmatic  significance  of  distinguishing  between  basic  versus  local 
planning.  It  is  known  (see  references  in  [Yap85])  that,  even  for  the  simplest  classes  of 
robots,  when  the  degrees  of  freedom  of  the  robots  is  arbitrarily  large  then  motion- 
planning  is  probably  intractable  (NP-hard  or  worst).  Even  if  the  number  of  degrees  of 
freedom  is  fixed  (e.g.  rigid  polygonal  robots),  the  complexity  of  motion-planning  of 
arbitrarily  complex  robots  seems  to  require  high  degree  polynomials;  this  complexity  is 


The  interested  reader  may  refer  to  the  Strang's  paper  for  the  amusing  genesis  of  the  problem. 


unacceptable  in  practice.  However,  a  practical  motion-planning  system  can  exploit  the 
heuristic  fact  that  'most'  robot  motion  takes  place  in  relatively  uncluttered  space.  In 
uncluttered  space,  the  robot  B  can  be  approximated  by  some  simpler  enclosing  shape  B' . 
For  example  if  we  enclose  B  in  the  minimal  size  disk  B'  ,  then  fast  O(filogn)  algorithms 
are  available.  At  the  isolated  localities  that  are  cluttered,  and  only  then,  we  invoke  more 
sophisticated  algorithms  that  take  the  detail  shape  of  B  into  account.  A  further  key 
observation  is  that  most  of  these  cluttered  localities  are  reasonably  stereo-typed,  such  as  a 
door  or  a  corner  in  a  corridor.  This  is  especially  true  in  a  controlled  or  planned 
environment  such  as  a  factory  floor.  We  should  be  able  to  provide  local  experts  for  such 
stereo-type  situations.  This  paper  is  therefore  one  conaibution  towards  providing  such 
local  experts. 

The  rest  of  this  paper  is  organized  as  follows:  section  2  gives  some  preliminary 
concepts.  In  sections  3  and  4  we  define  three  notions  of  width.  We  call  these  the  door- 
width,  ring-width  and  elastic  ring-width,  respectively.  This  is  partly  to  indicate  that  for 
non-convex  polygons,  the  concept  of  width  is  not  unique  and  can  be  subtle.  The  main 
technical  results  of  the  paper  are  in  sections  5,  6  and  7  where  we  develop  an  algorithm 
for  computing  the  door-width  of  a  polygon.  We  conclude  in  section  8.  There  is  an 
appendix  containing  some  geometric  lemmas. 

2.    Preliminaries 

In  this  paper,  we  want  to  refer  to  a  'physical'  chair  (a  polygonal  body)  moving  in 
the  plane:  by  'physical'  we  mean  that  the  chair  is  a  geometric  figure  that  has  no  fixed  or 
instantiated  position  in  the  plane.  I'o  make  this  idea  precise,  we  must  distinguish 
between  two  coordinate  systems  (or  'frames  of  reference'):  the  standard  frame  and  the 
actual  frame.  Points  in  standard  (resp,  actual)  frame  will  be  designated  as  standard  (resp. 
actual)  points. 

A  chair  (or  bodv)  fl  is  a  compact  set  of  standard  points  whose  boundary  S{B)  is  a 
simple  closed  polygonal  path.  Throughout  this  paper,  we  fix  some  arbitrary  B  and  let  n 
be  the  number  of  corners  of  B.  Typically  but  not  necessarily,  the  origin  of  the  standard 
frame  is  a  point  in  B.  .Note  that  in  general,  if  we  treat  several  bodies  then  each  body 
must  be  given  its  own  standard  frame. 

It  IS  important  to  note  that  there  are  only  conceptual  but  no  formal  distinctions 
between  the  standard  and  actual  points  However,  we  do  not  mix  concepts  from  the  r^vo 
frames.  For  instance  if  x  is  a  standard  point  and  K  is  a  set  of  actual  points,  the 
expression  'x  €  Y'  is  not  meaningful  Or  again,  if  X  is  a  set  of  standard  points,  we  would 
not  use  the  expression  '.YPiK'  Instead,  the  connection  between  the  two  frames  is  always 
done  via  certain  transformations,  as  we  now  explain:  A  placement  (of  the  standard  frame 
m  the  actual  frame)  is  represented  bv  a  rigid  transformation  of  the  standard  frame, 
specified  by  a  triple  Z  ^  (.r,  y,  0)  where  (x .  v)  is  a  pair  of  reals  and  0  <  6  <  2t7.  Ihc 
placement  Z  specifies  that  the  origin  of  the  standard  frame  is  at  the  actual   point  (x ,  y) 


and  the  positive  x-axis  of  the  standard  system  is  to  be  given  the  orientation  6.-  More 
precisely,  for  any  standard  point  p,  the  transformation  Z  =  (x,  y,  6)  of  p  is  the  actual 
point  defined  by 

p[Z]  =  (jc +  p  jCose -p^sine,  y  *p,sine  -  p,cos6) 

which,  in  matrix  notation,  is 

where  cis{Q)  is  a  2x2  rotation  matrix.  If  5  is  a  set  of  standard  points  then 
S[Z\  =  {p[Z]  :  p  €  S) .  For  instance  B[Z]  denotes  the  set  of  actual  points  occupied  by  B 
at  placement  Z.    Also  let  9(Z)  denote  the  orientation  of  Z. 

The  [Z]  and  <Z>  notations:  For  each  Z,  we  can  regard  [Z]  as  a  map  from  the 
standard  frame  to  the  actual  frame  of  reference  It  is  convenient  to  have  a  notation  for 
the  inverse  map:  if  p  is  an  actual  point,  define  p<Z>  to  be  the  standard  point  q  where 
q[Z]  =  p.    Thus  <Z>  is  the  inverse  of  the  map  [Z\. 

The  set  D?  of  all  placements  is  essentially  £-x  5'  where  £-  is  the  Euclidean  plane  and 
5'  the  unit  circle.  (D?  is  also  called  the  Euclidean  group  or  rigid  planar  transformations.) 
We  can  make  Df  into  a  metric  space  K(fi)  by  defining  the  distance  d{Z,  Z')  between  two 
placements  Z ,  Z'  to  be  the  Hausdorff  distance  between  B[Z]  and  B[Z'\.  Note  the 
dependence  of  this  metric  on  B. 

Let  W  be  a  set  of  actual  points  called  the  obstacle  set.  We  define  Z  to  be  free  (for  B  , 
with  respect  to  W)  if 

interwr{B[Z\)r\W   =   0. 

(Note:  in  the  literature,  the  alternative  condition  fl[Z]nW  =  0  is  sometimes  used  to 
define  free  placements.)  The  free  space  of  B  with  respect  to  H'  is  the  set  of  2  that  are 
free;  let  FP(W,  B)  C  iR(B)  denote  this  free  space.  If  B  is  understood,  we  write  FP{W) 
instead  of  FP{W ,  B).  A  motion  of  B  amidst  obstacle  \V  is  a  continuous  path 
IT  :  [0,  1]  -  '}?(fl).  The  motion  is  free  or  (avoids  the  obstacle  W)  if  tt(;)  €  FP(W,  B)  for 
all  t.  We  sometimes  denote  ~(t)  by  Z,  and  write  -  =  {Z,  :  r  t  [0,  1]}  C  J?(fl).  A 
constant  motion  tt  is  one  in  which  7t(0  =   -(0)  for  all  t. 

Given  t^A'o  motions  -^^  and  tt,  where  t7,-(1)  =  t7,(0),  we  define  their  composition  :t 
by: 


7T(r)  = 


-iTo(20         if  t  <  1/2 


IT. (2;-  1)    if  (  >  1/2 


•  .\  vivid  picrure  comes  irom  thinking  of  ;he  actual  trame  as  the  table  on  which  we  place  a  sheet  ot 
paper  (the  standard  frame)  that  has  a  drawing  ot  the  chair. 


We  denote  this  composition  by  tt  =  tto'.tti. 

Two  motions  ttq  and  -n ^  are  said  to  be  similar  if  there  is  a  continuous  monotonic 
function  /  :  [0,  1]  -  [0,  1]  such  that  /(O)  =  0,/(l)  =  1  and  TTgtf)  =  -iT|(/(r))  for  all  r. 
Otherwise,  ttq  and  tt,  are  dissimilar.  For  instance,  it  is  not  hard  to  see  that  -n Q;{-n ^ ,t! ^) 
and  (iTQi-iT,);-!!-,  are  similar  provided  the  compositions  in  these  expressions  are  well- 
defined.  We  distinguish  motions  only  up  to  similarity,  and  hence  we  will  be  happy  not  to 
parenthesize  a  sequence  of  compositions:  the  above  will  be  written  as  ttq-.v ^-.-n^- 

3.    The  Door  Width 

We  formalize  the  first  notion  of  the  width  of  a  body  B .  Let  w  >  0  and  consider  the 
following  set  of  actual  points: 

0„   =   ((0,  >)  :  y  <  0  or  v  >  w}  C  £- 

We  call  D,^  a  door  of  width  w.  This  captures  the  notion  of  an  infinite  wall  with  a  hole  of 
width  w.  It  is  convenient  to  have  a  notation  for  the  following  partition  of  the  actual 
plane  by  the  v-axis: 

LHS  =   {{x,  y)  -.x  <  0}  ,       ('left  hand  side') 

RHS  =   {(jt,  y)  -.x  >  0}  ,       {'right  hand  side') 

Y  =  {{x,y)  ■.x  =  0),       i'y-axis') 

We  say  a  motion  (Z,  :  f  €   [0,  1]}  of  fl  is  a  passage  through  D     if  it  avoids  D  ^  and  also 

B[Z^]  C  LHS     and     B[Z,]  C  RHS. 

In  other  words,  the  body  moves  from  the  left  half-plane  to  the  right  half  while  avoiding 
collision  with  the  wall. 

Definition.    The  door-width  of  B   is  the  smallest  w  such  that  there  is  a  passage  through 

This  is  equivalent  to  the  definition  of  Strang  in  case  B  is  convex.  The  door-width  of 
a  convex  body  is  simple  to  compute:  let  (v,  e)  denote  a  (vertex,  edge)  pair  of  R  .  We  say 
(v,  e)  IS  a  opposite  pair  if  the  line  through  p  and  parallel  to  f  is  a  supporting  line  ol  B. 

Lemma  1. 

(a)  The  door-width  of  B   is  attained  at  some  opposite  pair  {v ,  e]  m  the  sense  that  the 
door-width  IS  equal  to  the  distance  between  v  and  the  line  through  e 

(b)  There  are  0(n)  opposite  pairs  and  these  can  be  obtained  m  linear  time 

Part  (b)  is  easy,  and  part  (a)  follows  from  Strang's  characterization  [StrSZ]  that  the 
width  of  any  convex  shape  S  is  obtained  as  the  projection  of  S  onto  some  line.  The 
natural  extension  of  this  characterization  to  3  dimensions  is  false,  and  this  is  one  of  the 


points  of  Strang's  paper.    The  algorithmic  implication  of  this  lemma  is  that  we  need  not 
rotate  a  convex  chair  inside  the  doorway. 

It  follows  that  the  door-width  of  convex  bodies  can  be  computed  in  linear  time. 
Now  consider  the  following  body  Sq  consisting  of  a  2x  (2^  e)  rectangle  with  a  Ix  t  slit  on 
one  of  the  sides  of  length  2  +  e.  Here  1  >  €  >  0  and  it  is  not  hard  to  see  that  Bq  has  a 
door-width  of  2: 

2 


V2 


/ 


/ 


2^€ 


Figure  1.  Body  fl(- 


4.    The  Ring  Width  and  Elastic  Width 

The  purpose  of  this  section  is  to  show  that  the  notion  of  width  for  non-convex 
objects  admits  other  reasonable  interpretations.  This  section  may  be  omitted  since  later 
sections  are  independent  of  it.  The  reader  may  intuitively  feel  that  the  width  of  Bq  in  the 
previous  example  could  be  defined  so  that  it  is  V'2  rather  than  2.  This  intuition  comes 
from  imagining  a  'ring'  of  diameter  \''2  that  is  passed  over  the  body  B,^.  However,  the 
correct  definition  of  what  it  means  for  a  body  to  pass  through  a  ring  is  a  little  subtle  (the 
reader  is  invited  to  attempt  such  a  definition).  For  instance,  a  suitable  definition  must 
contend  with  the  following  example;  a  ring  of  width  1  can  superficially  sweep  every  point 
of  the  following  non-convex  quadrilateral  even  though  the  ring  cannot  pass  through  the 
body: 


Figure  2.  Non-convex  quadrilateral 


To  capture  the  concept  precisely,  first  define  the  ring  of  diameter  w  >  0  to  be  the 
pair  of  actual  points 

/?„,  =  {(0,0),  (0,w)}. 

It  is  also  convenient  to  have  a  notation  for  the  (actual)  line  segment 

A.  =   {(0,  >•)  :  0  s  V  s  w}  C  £-. 

Consider  a  motion  -rr  =  {Z,  :  r  €  [0,1]}  of  B  that  avoids  /?„  such  that  B[Z^]  C  LHS  and 
B[Z^]  C  RHS.  For  any  standard  point  p  i  interior(B),  let  J^  C  [0,  1]  be  the  set  of  time 
instances  i  such  that  p[Z,]  lies  in  /„  .  Note  that  well-known  properties  of  the  unit  interval 
imply  that  J    is  a  finite  union  of  disjoint  intervals  of  the  form 

Jp  =   [^1-  'i]U[i:,  ':]U U[5,,;,J,     (m  >  1) 

where 

0  <  ii  ^  r,  <  i,  <  f,  <    ■  •  •    <  j^  <  r„  <  1. 

If  t  is  one  of  the  f,'s,  and  s  is  one  of  the  i,'s,  define 

■l;2  ifp[Z,_]CLW5 


X~(^)  =   X'is,-) 


1/2  othenx'ise 


x"(0  =  x"('.  ^) 


1/2  ifpfZ,.]  C  /?W5 

-  1/2  otherwise 


The  notation  Z,_  in  the  preceding  definition  refers  to  the  placement  of  B  at  some  moment 
preceding,    but   'sufficiently    close'    to,    /.     Similarly    for   Z,..     Now    we    define    the    -n- 


crossing  number  of  p  €  interior{B)  (relative  to  it)  to  be 

x(p)  =  x(p-^)  =  ix-(^,)  -  Ix'iO- 

.-1  1=1 

Essentially,  the  crossing  number  of  p  increments  (resp.  decrements)  by  one  each  time  p 
crosses  from  LHS  (resp.  RHS)  to  the  RHS   (resp.  LHS).     It  is  easy  to  see   that  since  B 
moves  from  LHS  to  RHS  that  the  crossing  number  of  any  interior  point  p   must  be  an 
integer  value. 
Definition.    A  motion  it  =  {Z,  :  r  €   [0,  1]}  is  a  passage  through  /?,^  if  it  avoids  /?,,,  and 

(1)  B[Zq]  C  LHS  and  fl[Z,]  C  RHS,  as  before. 

(2)  For  each  point  p  in  interior{B),  the  u-crossiiig  number  of  p  is    >  1. 

It  can  be  shown  that  for  any  passage  tt,  the  Ti-crossing  number  does  not  depend  on 
the  choice  of  the  interior  point  p  (use  the  fact  that  B  is  connected). 

Definition.    The  ring-width  of  B  is  the  smallest  w  such  that  there  is  a  passage  through  R^. 

It  is  now  easy  to  see  that  the  ring-width  of  Aq  in  the  last  example  is  indeed  Vl. 

There  is  one  question  which,  though  obvious  for  door-width,  is  not  clear  for  ring- 
width.  Does  the  ring-widths  satisfy  the  following  monotonicity  property:  if  there  is  a 
passage  of  B  through  a  ring  of  diameter  w  then  there  is  a  passage  of  B  through  any  ring 
of  diameter  w'  where  w'  >  w. 

Conjecture  1.    Ring-widths  have  the  monotonicity  property. 

Another  possible  notion  of  width  is  what  we  call  an  elastic  width.  This  is  essentially 
the  same  as  the  ring-width  except  that  we  have  an  elastic  ring  that  (can  continuously) 
shrink  to  a  smaller  width  if  necessary.  We  omit  a  formal  definition.  It  is  clear  that  for 
any  body  B , 

door-width(B)  >  ring-width(B)  2  elastic-width(B) . 

It  is  clear  that  these  two  inequalities  are  in  fact  equalities  for  a  convex  B .  Our  example  Sq 
above  shows  that  the  first  inequality  may  be  strict.  We  conjecture  that  elasticity  does  not 
help: 

Conjecture  2.    For  any  chair  B  ,  rifi^-uit/r/!  (S)   -   elastic-width{B). 

5.    Constraint  analysis 

Recall  that  Y  is  the  v-axis  of  the  actual  plane.  For  any  Z  i  iR(fl).  if 
inierior(B[Z])r\Y  7^  0  we  say  that  Z  has  positive  clearance,  otherwise  it  has 
zero  clearance.  If  Z  has  positive  clearance,  there  is  a  unique  pair  of  standard  points 
(u,  v)  such 

interioriB[Z])nY  C  {u[Z],  v[Z]) 


and  the  distance  between  u  and  v  is  minimum.    Assuming  that  u[Z\  lies  above  v[Z]  on  Y , 
we  use  the  notation 

top<Z>         and         boi<Z> 

to   refer  to   these   points   u    and   v,   respectively.     If  Z   has  zero  clearance,   top<Z>    and 
boi<Z>  are  undefined. 

Y 


Figure  3.  Illustrating  the  definition  of  clearance 


For  any  placement  Z,  define  Clearance{Z)  to  be  0  if  Z  has  zero  clearance;  otherwise 
Clearance(Z)  is  the  distance  between  top  <Z>  and  bot<Z> .  For  any  motion  -,  let 
Clearance{-n)  be  the  maximum  of  Clearance (-rrii)),  t  €   [0,  1]. 

Two  placements  Z  and  Z'  are  equivalent  if  they  are  vertically  translatable  to  each 
other:  more  precisely,  if  Z  =  (x,  v,  6)  and  Z'  =  (j:',y'.9)  then  r  =  x'  and  6  =  6'. 
Two  motions  -  and  tt  '  are  equivalent  if  some  motion  tt''  similar  to  tt  is  pointwise 
equivalent  to  tt".  i.e.,  for  all  t,  7r"(f)  and  ■tt'(;)  are  equivalent.  Despite  the  asymmetry 
in  this  definition,  the  relation  is  a  true  mathematical  equivalence  relation.  Note  that  two 
equivalent  placements  or  motions  have  the  same  clearance.  Again,  we  normally 
distinguish  motions  up  to  equivalence. 

Definition.     A   pair   (e ,  e' )  of  distinct  edges  of  5    is  called  a  bridge   if   there   exists   a 
placement  Z  with  positive  clearance  such  that 

(a)      Y<Z>  contains  no  vertices  of  B  and 


(b)     top<Z>  and  boi<Z>  (resp.)  lie  in  e  and  e' . 

We  say  such  a  placement  Z  is  on  the  bridge  (e,  e').  Let  5  C  9J(fl)  be  the  set  of  all 
placements  on  the  bridge  (e,  e').  The  closure  of  S  in  X{B)  is  called  the  channel  of  {e ,  e') 
and  denoted  C(f,  ?'). 

Observe  that  because  of  (b)  in  the  definition  of  a  bridge,  we  regard  (e,  e')  and 
(e',  e)  as  distinct  bridges.  Indeed,  bridges  and  hence  channels  always  come  in  such  pairs. 
Let  V(e,  e')  denote  the  set  of  vertices  of  B  that  lies  outside  the  interior  of  the  rectangle 
(or  a  triangle  in  the  degenerate  case)  formed  by  the  endpoints  of  e  and  e' .  Thus  the  3  or 
4  distinct  endpoints  of  e  and  e'  are  included  in  V(e ,  e').  U  Z  is  on  the  bridge  {e ,  e')  then 
the  y-axis  Y<Z>  partitions  V(e,  e')  into  two  non-empty  sets  L  and  R.  It  is  not  hard  to 
see  that  iL,R)  and  ie,e')  each  uniquely  the  other.  Hence  we  sometimes  denote  the 
channel  C(e ,  e')  by  C(R ,  L).  Further  observe  that  a  channel  forms  a  connected  subset  of 
5}(B). 

Form  the  respective  convex  hulls  H^  and  W^  of  L  and  R.  A  point  p  on  the 
(boundary  of  the)  convex  hull  H^  or  //;,  is  said  to  be  visible  in  the  channel  C{L,  R)  if 
there  is  a  placement  Z  in  the  channel  such  thatp[Z]  lies  on  the  y-axis,  i.e.,  p  i  Y<Z>. 
We  conventionally  arrange  for  edges  e ,  e'  of  a  bridge  to  have  an  east-west  orientation 
with  H [^  on  the  west  and  H p,  on  the  east.  The  set  of  visible  points  of  //^  (resp.  H^)  form 
a  contiguous  chain  on  (the  boundary  of)  H^  (resp.  //^)  which  we  call  the  west  bank  (resp. 
east  bank).  A  placements  Z  €  C(f ,  <?')  is  said  to  be  on  the  east  bank  if  Y <Z>  intersects 
the  east  bank.  Similarly  for  the  west  bank.  It  is  easy  to  see  that  the  boundary  of  a 
channel  (as  a  subset  of  R{B))  is  composed  of  precisely  the  placements  on  the  two  banks. 
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Y<Z> 

Figure  4.  Vertex  v  lies  on  the  east  banlc  of  the  {e ,  e')  channel 

Our  main  problem  in  this  section  is  show  how  to  obtain  canonical  motions  in  channels. 

Definition.    Let  -u  =  {Z^  :  t  i   [0,1]}  be  a  motion  of  B  .    We  say  tt  is  a  pure  motion  if  it  has 
one  of  these  two  forms: 

(i)  TT  is  a  pure  rotation  if  there  exists  a  standard  point  p  such  that  for  all  t,  p[Z]  is 
constant  and  6(Z,)  is  monotonic  (increasing  or  decreasing)  in  r.  The  point  p  is 
called  the  center  of  the  rotation. 

(ii)  TT  is  a  pure  translation  if  6(Z,)  is  a  constant  for  all  t  and  for  each  standard  point  p, 
p[Z]  moves  raonotonically  in  a  straight  line.- 

We  would  like  to  further  reduce  pure  motions  in  a  channel  to  a  very  simple  form. 
Translations  have  simple  geometry  so  the  main  work  lies  in  studying  pure  rotations  in  a 
channel.  To  do  this,  consider  the  following  setup:  fix  a  channel  C(e.  e' )  and  let  tt  be  a 
pure   rotation   in  the  channel  centered   at  a  visible    vertex    v   of  the   channel      We   are 


rigid. 


'  It  IS  enough  that  any  r»o  standard  points  p    ^  p'  move  in  •>traight  lines,  since  transformations  are 
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interested  in  the  clearance  function  £>(/)   =   Clearance (iric)).    In  the  appendix,  we  show 

Lemma  2.    The  function  D(l)  has  at  most  one  maxima  t   =   tr-^  (:   (0,   1) 

If  f,,  exists,  wc  say  that  the  orientation  of  the   line   Y  <t^(tf^)>    is  special  for  v   relative  to 

To    define    the    pure    motions    of    interest    we    introduce    the    auxiliary    concept    of 
'constraints'. 

Definition.  Let  C{e,  e' )  be  a  channel.  A  constraint  A^  is  a  syntactic  (i.e.  formal)  object 
associated  with  a  channel  C{e,  e')  with  one  of  the  following  four  forms: 

K{e),  K{e'),  K(v),  K  {u ,  v),  A'(v,  e.  e' ) 

where  u  ,v  are  visible  vertices  of  the  channel  C(e,  e').    We  say: 

(A)  Z  satisfies  K{v)  if  v[Z]  €   Y  and  v  is  a  visible  vertex  in  C  {e ,  e'). 

(B)  Z  satisfies  K{e)  if  the  y-axis  Y  is  perpendicular  to  c[Zj.    Similarly  for  K(e'). 

(C)  Z    satisfies   K{u,  v)    if   Y   is    parallel    to    the    line    u[Z]v[Z]    where    u ,  v    are    adjacent 
vertices  on  one  bank  of  the  channel 

(D)  Z   satisfies   K{v,e,e')   if   the   orientation  of   Y<Z>    is  the  special   orientation  of  v 
relative  to  the  bridge  {e ,  e'). 

Each  constraint  K  is  associated  with  a  predicate  (also  denoted  K)  on  the  channel  where 
the  predicate  K  holds  at  Z  if  Z  satisfies  K  according  to  the  preceding  definition.  We 
often  confuse  the  distinction  between  the  constraint  and  its  predicate. 
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Y<Z> 


Y<Z> 


(A)    K{.) 


Y<Z> 


(C)    K(u,v) 


(D)    K(v^e,  €■) 


Figure  5.  Four  types  of  constraints 


A  constraint  K  is  trivial  if  the  corresponding  predicate  is  identically  false.  (Clearly 
no  constraint  can  be  identically  true.)  For  instance,  suppose  the  number  of  visible  vertices 
of  a  channel  is  k.  Note  that  k  ^  .1.  Then  there  are  k  nontnvial  constraints  of  type  (A), 
at  most  2  nontnvial  constraints  of  tvpe  (  B).  k  -  2  mm  trivial  constraints  of  type  ( C),  and  at 
most  k  nontnvial  constraints  of  type  (  D). 

Two  constraints  K,  .ind  A',  of  .i  channel  are  dependent  if  thev  arc  identical  as 
predicates,  otherwise  thev  are  tndepi-ndent  -     The  next  lemma  indicates  that  cimstraints  of 
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satisfy  some  other  constraint  K  of  the  form  K{u' ,  v')  or  K{e)  or  K{u,  e,  «').  We  will 
charge  the  critical  placement  Z  to  K .  Note  that  K{u' ,  v')  can  be  charged  0(1)  times,  and 
there  are  O(n')  constraints  of  the  form  K{u' ,  v').  The  lemma  holds  for  this  case.  Next, 
each  K{e)  is  charged  0{n)  times.  To  see  this,  each  e  is  involved  with  at  most  n  other 
edges  to  form  channels.  But  each  channel  causes  K{e)  to  be  charged  0(1)  times.  Since 
there  are  overall  n  such  choices  for  K(e),  the  result  also  follows.  Finally,  each 
K{u,e,e')  is  charged  0(1)  times;  it  remains  to  show  that  there  are  O(n-)  such 
constraints.  Each  constraint  of  the  form  K(u,  e,  e')  can  be  charged  to  u;  it  is  easy  to  see 
that  u  can  be  charged  only  0(n)  times  and  there  are  0{n)  choices  for  u.  The  proof  is 
now  complete.  Q.E.D. 
Definition.    A  pure  motion  -77=   (Z,  :  r  €  [0,  1]}  is  semi-elementary  if  the  following  holds: 

(i)  If  -i;  is  a  translation  then  there  exists  a  constraint  K  of  the  form  K{e)  or  K{u,  v)  or 
K{u,  e,  e')  such  that  for  all  /,  Z,  satisfies  K.  If  tt  is  a  rotation  then  there  exists  a 
constraint  of  the  form  K{v)  such  that  for  all  r,  Z,  satisfies  K . 

(ii)     Z,  is  not  critical  for  all  0  <  r  <  1. 

If,  in  addition,  Zg  and  Z,  are  critical  then  we  say  that  it  is  elementary. 

One  special  type  of  elementary  translation  tt  deserves  special  mention:  if  T7(f)  is 
equivalent  to  tt(0)  for  all  t,  we  call  it  a  slide. 

It  is  not  hard  to  see  that  for  any  two  elementary  motions  if  their  initial  and  final 
placements  are  respectively  equivalent  then  they  are  equivalent  motions.  In  other  words, 
up  to  equivalence,  elementary  motions  may  be  identified  with  the  equivalence  classes  of 
their  initial  and  final  placements.  This  shows  an  upper  bound  of  0(n"')  elementary 
motions.    We  improve  this  estimate  considerably. 

Theorem  6.    There  are  0{n^)  inequivalent  elementary  motions. 

Proof.  It  is  easy  to  see  that  each  critical  placement  Z  is  the  starting  placement  for 
at  most  3  inequivalent  elementary  translations,  namely,  a  slide  and  at  most  two 
translations  (in  the  2  directions  normal  to  y-axis).  The  final  placements  of  such 
elementary  motions  are  completely  determined.  Thus  the  0{n-)  bound  holds  in  the  case 
of  elementary  translations. 

To  bound  the  number  of  elementary  rotations,  we  consider  two  cases  of  critical 
placements  Z.  If  Y<Z>  contains  only  one  vertex  of  B  (there  must  be  at  least  one  vertex) 
then  it  is  easy  to  see  that  there  are  at  most  2  elementary  rotations  starting  at  Z.  Each 
such  rotation  can  be  charged  to  Z.  Each  such  Z  is  charged  at  most  twice  and  there  are 
O(n-)   choices   for  such  Z.     The  second  case   is  where   Y<Z>    contains   more   than  one 

vertex  of  B .    Suppose  Y <Z>  contains  the  vertices  u, Uj.,  for  some  A  2:  2,  occurring 

in   the  order   indicated.     Each   elementary    rotation   centered   at   m,    (/   =    1 k)   and 

starting  from  such  a  Z  can  be  charged  to  a  consecutive  pair  of  vertices,  either  (u,,  «,^ ,)  or 
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(m,_i,  m,),  arbitrarily.  Clearly  no  pair  (u,  v)  of  vertices  is  charged  more  than  4  times  and 
there  are  only  O(n^)  such  pairs.  Hence  there  are  O(n-)  inequivalent  elementary  rotations 
overall.    Q.E.D. 

6.    How  to  cross  a  channel 

A  nice  motion  -n  is  either  a  semi-elementary  motion  or  else  that  can  be  written  in 
the  form 

■n  =  ■T(o^,■n^■,  •  •  •  lir^       {k  ^  2) 

where 

(i)      ttq  and  tt^  are  semi-elementary  rotations, 

(ii)     iT|,  .  .  .  .TTt-i  are  elementary  motions  (translations  or  rotations). 

(iii)    No  other  motions  from  -it(O)  to  tt(1)  achieves  a  smaller  clearance. 

For  example,  elementary  motions  are  nice.  Note  that  -q(1)  and  -n ^(0)  are 
necessarily  critical. 

We  show  that  any  motion  can  be  reduced  to  a  nice  one.  We  begin  by  considering 
motion  in  a  fixed  channel  C{e,  e').  As  usual,  we  talk  about  the  east  and  west  banks  of 
the  channel.    We  define  a  map 

Im  :  C{e,  e')  -  C  (e ,  e") 

that  takes  any  placement  Z  €  C(e,  ?')  by  a  direct  translation  to  a  placement  Im{Z)  on  the 
west  bank.  .More  precisely,  for  any  Z,  let  v  be  the  vertex  on  the  west  bank  that  is  closest 
to  Y<Z>  If  d  is  the  distance  of  v  from  Y<Z>  and  Z  =  (.t,  v,  6)  then 
lm{Z)  =  {x~d,  V,  6).    In  particular,  v  d   Y<Im{Z)>. 
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west 


We  may  assume  that  the  bridge  is  narrower  on  the  west  bank  than  the  east  bank,  i.e., 
Clearance{Z)  2:  Clearance{Im{Z)) .  Given  any  motion  tt  in  a  channel,  we  define  a  motion 
/m(iT)  on  the  west  bank  io  the  natural  way. 

Lemma  7.    Let  Zq  be  on  the  east  bank  of  a  channel.    Then  there  is  a  nice  motion  r,  from 
Zq  to  Imi^Zo)  that  achieves  its  maximum  clearance  at  Zq. 

Proof.     We  construct  it  =  TTpi-n  pTT,  as  follows: 

TTj:     Rotate  Z^  in  the  direction  of  decreasing  clearance  until  the  first  critical  placement 
Z,.  (If  Zq  is  critical,  then  ttq  is  a  constant  motion  and  Zq  =  Z,.) 

■iT|:     Translate  from  Z,   across  the  channel  to  Im{Z^).    This  motion  is  elementary.     Note 
that  tT|.|;-:T|  has  non-increasing  clearance. 

ttt:     Finally,  rotate  from  Im(Z^)  directly  to  /m(Zo).    This  motion  is  semi-elementary  and 
hence  achieves  its  maximum  clearance  at  /m(Z,)  or  at  lm{ZQ). 

Since  the  clearance  at  lm{Z,)  and  ^m(Z-,)  are  both  no  larger  than  the  clearance  at  Z^,  we 
conclude  that  tt  attains  its  maximum  clearance  at  Zg.    Q.E.D. 
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Lemma  8.    Let  Zq  and  Z,  be  on  the  banks  of  a  channel.    Then  there  is  a  nice  motion  it 
from  Zq  to  Z I . 

Proof.      There  is  an  obvious   nice   motion  it  ,   from  lm{Z^)   to  Im(Z^)   consisting  of 
rotations  only.    Furthermore,  no  motion  from  Z,^  to  Z,  has  clearance  less  than 

Tmx{Clearance(ZQ),  Clearance{Z^),  Clearance{-u  ^)] 

The  previous  lemma  gives  us  nice  motions  ttq  and  ir^  from  Zq  to  ImiZ^)  and  from  Z,  to 
/m(Z,)  (respectively).  Since  the  clearance  of  Tj  and  -r,  are  attained  at  Zg  and  Z, 
(respectively)  it  follows  that  ttq;it,;tt2  achieves  the  minimum  clearance  of  any  motion 
from  Zq  to  Z,.  The  two  consecutive  semi-elementary  rotations  in  ttq;-?;,  can  be  replaced 
by  an  elementary  rotation  so  that  the  result  is  a  nice  motion  ttj.  Repeating  this  for 
TTjiiT;,  we  get  the  desired  result.    Q.E.D. 

We  are  now  ready  for  our  main  result: 

Theorem  9.  Let  Z^^  i  LHS .  Z^  i  RHS.  If  there  is  a  passage  of  B  from  Zq  to  Zi  through 
a  door  D^  then  there  is  one  of  the  form. 

TTQJTTpTTT 

where  tTq  and  tt;  have  zero  clearance  and  are  each  the  composition  of  at  most  three 
pure  motions,  it  ,  is  a  composition  of  elementary  motions  that  achieves  the  minimum 
clearance  (i.e.  the  door-widih  of  Bj. 

Proof.  Suppose  TT  is  a  passage  of  B  from  Zg  to  Z,  whose  clearance  is  minimum. 
We  may  suppose  that  there  is  an  open  interval  /  C  [0,  1]  such  that  Ti(r)  has  zero 
clearance  precisely  when  t  ^  I.  We  may  further  suppose  that  for  some  A  s  1,  there  are 
moments 

^0  ^  'o  <  ^1  -  'i  <    •  ■  ■    <  ^*  ^  '^ 

such  that  /  =  (ig,  f.)  and  for  all  t  i  [0,  1|,  7t(0  satisfy  a  constraint  of  type  (A)  iff 
f  €  [5,,  r,]  for  some  1  =  0,  .      .  .  *. 

These  s'%  and  r/s  divide  -  into  portions  which  we  now  analyze  in  turn.    First,  we 

know  by  the  previous  lemma  that  each  portion  of  the  form  -[;,,  s  _,](/--  0 *-l) 

can  be  replaced  by  a  nice  motion.  .Next  consider  ttI,?,:]  for  some  (i,  r)  =^  (5.;), 
'  =  0,  .  .  .  ,*.  It  is  not  hard  to  see  that  TT[i,  f|  can  also  be  replaced  nice  motion.  Using 
the  same  trick  as  before  of  replacing  two  consecutive  semi-elementary  rotations  by  an 
elementary  rotation,  we  get  a  nice  motion  tt-.  from  71,(0)  =  •!T(.Sp)  to  71,(1)  =   7T(r). 

Since  tt,  is  nice,  wc  may  write  is  as  7r^;7T«;7r^  where  rr<  is  a  composition  of 
elementary  motions,  tt^  and  tt^  are  semi-elementary  rotations,  and  Tr^(l)  and  77^(0)  arc 
critical.    It  is  now  easy  to  give  a  motion  tt^  (resp.  tt.)  composed  of  at  most  three  pure 
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motions  (translation,   rotation  and  translation)   from  Zq  to  ti^(I)-    (rcsp.   HjCO)  to  Z,). 
With  TT,  set  to  TTj,  we  see  that  TTg;^^;'^-,  satisfies  the  theorem.    Q.E.D. 

7.    Algorithmic  details 

The  basic  idea  is  to  construct  a  weighted  graph  whose  edges  correspond  to 
elementary  motions.  We  then  run  a  shortest  path  algorithm  (Djikstra's)  on  the  graph  to 
discover  a  minimum  clearance  path. 

For  any  simple  polygon  B ,  let  G(fl)  be  the  node-weighted  (undirected)  graph  where 
each  node  represent  an  equivalence  clasi  of  critical  placements  and  whose  edges  (or 
arcsY  correspond  to  elementary  motions  that  connect  adjacent  critical  placements.  (The 
slides  would  be  represented  by  self-loops  but  we  may  ignore  explicit  consideration  of 
them  here.)  The  weight  attached  to  a  node  is  equal  to  the  clearance  of  the  critical 
placement  represented  by  the  node.  Call  a  node  initial  (resp.  final)  if  its  critical 
placements  Z  satisfy  B[Z]  C  LHS  (resp.  B[Z]  C  RHS).  In  particular,  initial  and  final 
nodes  have  zero  weights.  We  further  introduce  into  G{B)  two  special  nodes  called  the 
source  and  the  sink  and  introduce  arcs  connecting  the  source  to  all  the  initial  placements 
and  connecting  the  sinic  to  all  the  final  placements.  The  weights  of  the  source  and  sink 
are  zero. 

The  number  of  nodes  and  number  of  arcs  in  G{B)  is  O(n-)  where  B  has  n  vertices. 
Given  any  path  specified  as  a  sequence  {u^,  .  .  .  ,  u^)  of  nodes  in  G{B),  let  the  weight  of 
the  path  be  the  maximum  of  the  weights  of  u,  {i  =  1,  .  .  .  ,k).  From  the  main  theorem 
in  the  previous  section,  the  width  of  the  chair  B  is  the  minimum  weight  path  connecting 
the  source  and  the  sink.  The  minimum  weight  path  can  be  found  by  Djikstra's  algorithm 
in  0{n-)  steps. 

It  remains  to  show  how  to  construct  G{B)  in  0{n')  steps.  For  each  vertex  v  of  fl  we 
rotate  a  ray  R  about  v,  keeping  track  of  the  sequence  of  edges  that  intersect  the  rotating 
line.  In  actual  fact,  it  is  sufficient  to  keep  track  of  the  first  and  last  edges  intersected  by 
R  (these  edges  may  be  undefined).  As  will  be  shown  below,  this  list  of  first  and  last 
edges  for  each  v  is  sufficient  for  us  to  construct  G{B).  The  full  circle  of  2tt  radians  is 
thereby  broken  up  into  basic  intervals  during  which  the  rotating  ray  R  does  not  pass 
through  any  other  vertices.  This  process  is  easily  accomplished  in  Oin\ogn)  per  vertex  or 
O(n-logn)  overall.  However,  using  the  results  in  [Wel85],  [AAGHI86],  this  bound  can 
be  improved  to  0{n~)  time.  For  completeness,  we  briefly  show  how  to  do  this  assuming 
that  we  have  computed  for  each  vertex  v  of  fl ,  the  (circular)  list  L{v)  of  all  the  other 
vertices  in  the  order  encountered  by  the  ray  R  sweeping  a  full  circle  about  v.  It  is  known 
[CL85],[AAGHI86]  that  the  n  lists  L(v)  can  be  obtained  in  O(n-)  time.  For  simplicity  in 
this  discussion,  we  assume  that  no  three  vertices  are  collinear. 


'  To  avoid  confusion  between  the    edges'  of  the  polygon  fl  and  the  'edges'  of  G{B).  we  call  the  laner 
arcs     Similarly, 'vertices'  refer  fl  and  'nodes'  refer  lo  G  (fl ). 
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We  first  break  up  L(v)  into  two  sublists  L^{v)  and  Lg(v)  (resp.)  consisting  of 
vertices  encountered  by  the  ray  R  above  and  below  (resp.)  the  horizontal  line  through  v. 
Let  0^(v)  and  0g(v)  be  the  corresponding  list  of  orientations  of  the  ray  R  through 
vertices  in  L^{v)  and  Lgiv),  respectively.  For  any  orientation  9,  we  let  pred  (6)  (resp. 
succ,(Q))  to  refer  to  the  largest  (resp.  smallest)  orientation  6'  €  0(v)  where  6'  <  9 
(resp.  9  <  9').  For  example,  each  'basic  interval'  defined  above  is  of  the  form 
(9,  jucc,(9))  where  9  €  0(v).  For  each  9  €  Q(v)  we  define  /,(9)  to  be  the  first  edge  in  B 
seen  by  the  ray  from  v  at  any  orientation  strictly  between  9  and  5ucc,  (9).  If  this  edge  is 
undefined  then  /,,(e)  =  ».  Let  ^^(v)  be  the  list  of  (/,(9)  :  9  €  0^(v)}  corresponding  to 
0,(v). 

Lemma  10.    The  lists  <t^(v)/<7r  all  vertices  v  can  be  computed  in  0{n-)  time. 

Proof.  Suppose  the  ray  from  v  to  another  vertex  u  is  at  orientation  9  €  0^(v). 
Then  it  is  not  hard  to  see/  (9)  is  equal  to  one  of  the  following: 

(a)  f,(pred,{Q)) 

(b)  One  of  the  edges  incident  to  u. 

(c)  f,ipred^{Q)). 

Furthermore,  it  is  easy  to  check  in  constant  time  which  is  the  case.  Let  v,,  v,,  .  .  .  ,  v,  be 
the  vertices  of  B  in  order  of  decreasing  y-coordinate.  We  will  compute  the  lists  0^(v,) 
for  (  =  1,  2,  .  .  .  ,n  (in  the  indicated  order).  Using  an  O(nlogn)  precomputation  (using 
a  plane-sweep  method),  it  is  seen  that  for  each  list  4'^(\)  we  can  determine  the  first 
value  in  the  list.  To  obtain  successive  entries  in  the  list,  the  above  cases  (a)-(c)  apply. 
Cases  (a)  and  (b)  are  easy;  in  case  (c),  it  is  clear  that  we  must  lookup  the  list  <t'^(v,)  for 
some  j  =   I,  .  .  .  ,i  -  I. 

We  claim  that  with  O(n')  preprocessing,  the  value  /  (pred  (9))  in  case  (c)  can  be 

obtained  in  constant  time.  To  see  this,  we  see  that  6  i  <i>  ^(v)  corresponds  to 
9  +  TT  €  <t'fl(v',).  Using  an  0{n-)  preprocessing,  we  may  assume  that  there  is  a  pointer 
link  from  the  occurence  of  9  in  0|(v)  to  the  occurence  of  9  -  tt  in  ^^(v,).  Next,  for 
each  9'  i  fl'g(v,),  we  may  assume  there  is  a  pointer  link  to  pred   (9--)  in  'I'^(v,);  such 

links  can  be  obtained  in  linear  time  by  a  merge-like  process  on  ^ ^{y  )  and  <t'^(w  )  The 
above  two  types  of  links  shows  that  case  (c)  can  also  be  handled  in  0(1)  time.  This 
concludes  our  demonstration  of  the  lemma.    Q.E.D. 

We  now  show  how  to  detect  all  the  bridges  (e,  e').  Let  u  and  v  be  adjacent  vertices 
on  one  of  the  banks  of  C{e,  e')  with  u  and  v  lying  on  opposite  sides  of  the  bridge  (see 
figure). 
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Figure  7.  Detecting  the  bridge  (e ,  e') 


Suppose  the  orientation  from  v  to  u  is  6.    Then  (up  to  interchanging  the  roles  of  u  and  v, 
etc)  the  four  following  conditions  hold: 

(i)  /,(6)  =  e.  Furthermore,  let  u'  be  the  point  such  that  the  orientation  from  v  to  u'  is 
succ^(B).  Then  all  of  the  edges  incident  to  v  lies  on  the  side  of  the  line  uv  opposite 
to  u' . 

(ii)     /.(9)  =  ==• 

(iii)  f,Jpred^{Q  ~  tt))  =  e' .  The  edges  incident  to  u  must  lie  on  the  same  side  of  the  line 
uv  as  the  edges  incident  to  v,  as  described  in  (i). 

(iv)    Aipred  (6--))   =  oc 

This  test  can  be  done  for  each  pair  (u,  v)  for  a  total  time  of  O(n-). 

For  each  bridge  {e,e'),  we  claim  that  the  above  process  may  be  modified  to 
discover  its  left  and  right  banks.  Referring  to  figure  7,  consider  two  adjacent  vertices 
V,  w  on  one  bank  of  C{e,  e')  with  v,  w  on  the  same  side  of  the  bridge.  It  is  easy  to  set 
up  conditions  on  <t(v)  and  <t>(w)  similar  to  (i)-(iv)  which  detect  that  v  and  w  are  adjacent 
vertices  on  one  bank  of  a  channel.  Hence,  by  linking  together  these  adjacent  pairs 
vertices  on  each  banks  of  a  channel,  we  get  a  representation  of  the  banks  as  adjacency 
lists.    We  leave  the  details  to  the  reader. 

Finally,  we  determine  all  the  arcs  of  G{B)  --  this  is  equivalent  to  identifying  all 
elementary  motions.  We  first  refine  the  basic  intervals  as  follows:  by  merging  the 
orientations  of  0^(v)  with  {--0  :  6  €  0fl(v)},  each  basic  interval  in  0^(v)  is  broken  up 
into  subbasic  intervals.  Similarly  we  can  break  up  the  basic  intervals  of  ©^(v)  into 
subbasic  intervals. 
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Each  subbasic  interval  /  of  v  corresponds  to  a  rotation  in  some  channel  C(e,  <?').  If 
V  is  a  visible  vertex  of  the  channel,  then  clearly  this  rotation  is  a  composition  of 
elementary  rotations  in  C{e,e').  We  show  how  to  further  break  up  such  /  into 
elementary  rotations.  It  is  easy  to  subdivide  /  into  smaller  intervals  at  those  orientations 
where  constraints  of  the  form  K{e),  K{e')  or  K{v,  e,  e')  are  satisfied.  It  is  slightly  more 
involved  to  break  up  /  at  orientations  where  it  satisfies  K{u,  «')  where  u,  u'  are  adjacent 
vertices  on  the  bank  of  (e,  e')  opposite  v.  To  do  this,  we  process  each  channel 
separately.  Take  the  two  banks  (represented  as  linked  lists)  of  a  channel  and  by  a 
process  similar  to  merging  rwo  sorted  lists,  we  can  obtain  the  orientations  in  each 
subbasic  interval  /  corresponding  to  constraints  of  the  form  K{u,  u).  This  completes  the 
subdivision  of  the  basic  intervals  into  elementary  rotations.  Furthermore,  the  elementary 
translations  across  that  channel  can  be  obtained  during  the  same  merge  process.  This 
concludes  our  proof  of  the  main  algorithmic  result: 

Theorem  11.  The  graph  G(B)  can  be  constructed  0{n-)  steps  and  the  minimum  weight 
path  can  be  found  in  0{n')  time. 

Corollary.  The  door-width  of  a  simple  n-gon  can  be  computed  in  O(n-)  time.  The 
same  algorithm  can  also  constructs  a  motion  of  the  polygon  through  any  door  of 
sufficient  width. 

Remark.  The  above  algorithm  may  be  iped  up  in  practice  as  follows.  A  point  p  on 
the  boundary  of  a  polygon  P  is  said  to  be  visible  if  there  is  a  ray  originating  at  p  that 
does  not  intersect  the  interior  of  P .  The  set  of  visible  points  of  P  is  composed  of  a  set  of 
polygonal  path-segments  where  each  path-segment  is  a  connected  subset  of  the  boundary 
of  P.  By  connecting  these  path-segments  in  the  natural  way,  we  obtain  what  mav  be 
called  the  visible  hull  of  P.  Clearly,  door-width(P)  =  door-width(VH{P))  where  \'H(P) 
denotes  the  visible  hull  of  P .  It  follows  easily  from  (LL84]  that  VH(P)  can  be  computed 
in  O(nlog^)  time.  Hence,  it  may  be  more  efficient  to  replace  the  polygon  by  its  visible 
hull  before  running  the  door  algorithm. 

8.    Conclusion 

(1)  Our  introduction  of  constraint  analysis  here  appears  to  be  a  general  and  powerful 
tool  in  studying  geometric  problems  of  this  nature.  It  is  intimately  related  to 
retraction  approach  of  motion  planning  [Yap85].  Planned  future  papers  will 
illustrate  other  uses  of  such  analysis. 

(2)  Can  our  0{n-)  complexity  for  door-width  be  improved  to  subquadratic  bounds? 
Our  algorithm  is  optimal  among  those  which  must  construct  the  graph  G{B). 
This  IS  because,  in  the  worst  case,  G{B)  has  size  a>(n-).  For  instance,  the  following 
polygon  has  amege(n-)  channels  and  elementary  motions:  take  a  regular  convex  n- 
gon  and  at  each  vertex  'add  an  arrowhead'  (see  figure). 
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Figure  8.  Arrowhead 

In  any  case,  subquadratic  bounds  may  be  possible  for  special  chairs  (eg.    rectilinear 
ones). 

(3)  It  would  be  nice  to  prove  the  two  conjectures  concerning  ring-widths  in  section  4. 
We  expect  some  of  the  ideas  of  the  present  paper  can  be  extended  to  give  an 
algorithm  for  computing  the  ring-width. 

(4)  Several  issues  are  raised  by  the  practical  approach  to  motion-planning  suggested  in 
our  introduction:  First,  how  does  one  partition  an  environment  E  into  parts  which 
are  'clutter-free'  and  parts  that  are  'cluttered',  assuming  some  suitable  definitions 
for  these  intuitive  terms.  We  refer  to  [SS85]  for  one  notion  of  clutter-free 
environments.  Next,  given  a  method  of  partitiiniing,  we  need  a  systematic 
cataloguing  (stereotyping)  of  the  tight  situations.  Finally,  for  those  common 
stereotypes,  we  need  efficient  local  algorithms.  Besides  the  door  stereotype  of  this 
paper,  another  interesting  situation  is  a  corner  in  a  corridor  or  hallway.  In  this 
context,  we  refer  to  several  interesting  articles  on  the  so-called  'moving  sofa' 
problem  cited  in  [Str82]. 
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APPENDIX:  Behavior  of  the  clearance  function  during  a  rotation 

Our  goal  is  to  prove  lemma  2  in  section  4.  Let  C{L,  R)  =  C{e,  e')  be  a  channel. 
We  investigate  the  clearance  during  a  pure  rotation  it  in  the  channel.  Let  the  vertex 
V  6  Z.U/?  be  the  center  of  the  rotation.  We  may  suppose  that  e  and  e'  are  not  parallel, 
since  otherwise  the  analysis  would  be  trivial.  Let  A,  A'  be  the  (standard)  rays 
originating  from  a  common  point  0  and  containing  e,  e'  respectively.  Thus  v  lies  outside 
the  wedge  AA'. 

Choose  O  to  be  the  origin,  A  to  be  the  positive  x-axis,  A'  to  lie  above  the  x-axis. 
Wlog,  V  lies  below  the  x-axis  and  to  the  right  of  A'.  Let  Q  (resp.  Q')  be  the  foot  of  the 
perpendicular  from  v  to  the  line  containing  A  (rpsp.  A').  For  any  angle  6,  let /?g  denote 
the  ray  emanating  from  v  at  an  angle  6  measured  counter-clockwise  from  the  ray  v  ~  Q' . 
In  particular,  Rq  denotes  the  ray  v  ^  Q'  and  R_^  denotes  the  ray  v  -  g.    See  the  figure. 


Figure  8.  Case  a  £  t7/2 


Let  0  denote  the  range  of  orientations  where  ^g,  6^0,  intersects  both  A  and  A'. 
Clearly  the  pure  rotation  it  may  be  parametrized  by  6  €  0.  Define  D(6)  to  be  the 
distance  between  the  two  points  of  intersection  of  R^  with  A  and  .\'.  Clearly  lemma  2  is 
a  consequence  of  the  following. 

Lemma.     The  function  0(9)  is  either  monotonic  or  else  it  has  a  unique  maxima  and  a 
unique  minima  in  its  domain  0  of  definition . 

Proof.     Let  a  >  0  be  the  angle  between  A  and  A'.    First  assume  that 

0  <  a  <  -IT/2 

as  in  the  above  figure.  Let  H  =  \Q'  -  v|  be  the  distance  from  v  to  line  containing  .\'  and 
h  =   IQ  -  v|  be  the  distance  of  v  below  the  x-axis.    Then  it  is  easy  to  see  that 

o(e)  = r-  (1) 

cose     cos(e-a)  ^  ' 

Clearly  D(e)  s  0  for  9  f  0.  Let  0  =  [9.,  9,].  By  breaking  the  analysis  into  three  parts 
if  necessary,  we  may  assume  that  one  of  the  following  three  cases  hold: 
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(i)      -77/2  <  e,  <  6;,  <   -a 
(ii)    -as;e,<e-,<0 
(iii)  0  <  e,  <  9,  <  ir/2-a. 

We  shall  use  primes  (e.g.  D')  to  denote  differentiation  with  respect  to  6.    We  have 


D'(e) 


H 


cose 


tanO- 


cos(e  +  a) 


tan(e  T  a). 


(2) 


In  cases   (i)   and  (ii),  we   show   that  D'(9)  <  0  so   that  D   is  strictly  monotonic:  (i)   If 

-■n/2  <  e  <   -a,     this     follows     from     the     fact    that  >   ~ >  0     (since 

cose  cos(e*a) 

D(e)  2:  0)     and     tanS  <  tan(e  +  a)  <  0.       (ii)     If     -a:<Qs,0,     then     we     see     that 
tan6  s  0  s  tan(6  +  a),  so  again  D'  is  negative. 


Case  (iii)  is  unfortunately  more  intricate.    Setting 

f(Q)  =    sec(e  -  a)tan(e  ^  a) 
sec9tan9 


we  may  write 


We  will  show  that 


D'(e)  =   sec9tane[//-;iF(e)]. 


F"(e)  >  0. 


(3) 

(4) 
(5) 


This  implies  that  F{Q)  has  at  most  two  zeroes.    The  expression  (4)  then  shows  that  D'{6) 
also  has  at  most  two  zeroes.    We  show  that  D'(9)  cannot  have  exactly  one  zero.    This 
follows  from   the  fact  that  D'(0)  <  0  (from   (2))  and  D'{e.)  <  0  (since,  by  definition 
o(e,)  =  0). 

We  can  now  see  the  global  behaviour  of  D(9):  D(e)  is  monotonically  decreasing  in 
the  range  (i)  and  (ii),  and  possibly  after  an  interval  of  increase  in  range  (iii),  it  finally 
decreases  to  zero.    Hence  the  lemma. 

To  show  (5),  we  take  logarithms  in  (3), 

logF  =  [logsec(9^  a)*  logtan(6-  a)]-  [logsecO-  logtane]. 

Differentiating: 

F' 

—  =   [tan(9  -  a)  -  tan9]  -  [sec(e  -  a)cosec(9  -  a)  -  sec9cosecel.       (6) 

Differentiating  again, 


=  [5ec-(e  +  a)-sec-9J^  [(sec-(9  -  a)  -  cosec-(9  -  a))  -  (sec-S  -  cosec^e)] 


F" 
F 


F 


2[sec-(9  -  a)-  sec-9]-  [cosec-9  -  cosec-(9  -  q)] 


(7) 


Since  F(Q)  >  0  (from  (3))  and  it  is  easy  to  sec   that  the  RHS  of  (7)   is  positive,   this 
concludes  our  proof  of  (5). 
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It  remains  to  consider  the  case  where 

a  >  -it/2. 


Figure  9.  case  a  >  -rr/Z 


The  expression  (1)  for  0(6)  still  holds  except  that  the  range  0  now  satisfies 

-  tt/2  <  9  <  0  <  (fc  +  a)  <  till 

for   8  €  0.     Hence   tan9  <  0  <  tan(e-^a)    and   we   conclude  D'(e)  <  0   as    in   case    (ii) 
above.    Q.E.D. 
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